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ABSTRACT 

The highly luminous (> 10'^'' erg s~^) supersoft X-ray sources (SSS) are beheved to be 
Eddington limited accreting white dwarfs undergoing surface hydrogen burning. The 
' current paradigm for SSS involves thermally unstable mass transfer from a 1 — 2Mq 

, companion. However this model has never been directly confirmed and yet is crucial 

for the evolution of cataclysmic variables in general, and for the establishment of SSS 
as progenitors of type la supernovae in particular. The key SSS is RX J0513. 9-6951 
which has recurrent X-ray outbursts every 100 — 200 d (lasting for ~ 40 d) during 
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Q ■ which the optical declines by 1 mag. We present the first XMM- Newton observations 

I of RX J0513. 9-6951 through one of its optical low states. Our results show that as 

C/3 . the optical low state progresses the temperature and the X-ray luminosity decrease, 

' behaviour that is anti-correlated with the optical and UV emission. We find that as the 

optical (and UV) intensity recover the radius implied by the spectral fits increases. 
The high resolution spectra show evidence of deep absorption features which vary 
during the optical low state. Our results are consistent with the predictions o f the 
^ . white dwarf photospheric contraction model proposed bv lSouthwell et all l)1996() . 

Key v^rords: binaries: close - stars: individual: RX J0513. 9-6951 - Magellanic Clouds 
- X-rays: stars. 
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1 INTRODUCTION 

The supersoft X-ray sources (SSS) are a class of X-ray 
objects discover ed by the Einstein X-ray Obser vatory in 
the late 1970s llLong. Helfand fc Grabelsky||l98ll) . During 
the 1990s the class of SSS has been c onsiderably en - 
larged through the 7? 05^ T All-Sky Survey jTrumpeill992l) . 
Archetypal SSS sources include GAL 83 and CAL 87, whose 
defining characteristic is their extremely luminous emission 
(up to Lboi ~ 10^* ergs"^) at very soft (< 0.5 keV) X-ray 
energies. As such, the SSS are rendered undetectable in the 
Galactic plane due to the high level of X-ray absorption, 
with known sources lying predominant ly in the Magellanic 
Glouds and M31 (see e.g. lGreineilll99dl . 

That the SSS are low mass X-ray binaries (LMXBs - e.g. 
Ivan Paradiis fc McClintock 1995) was suggested soon after 
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their discovery by their optical spectra. Among the ROSAT 
discoveries are several systems which also show these hall- 
marks, incl uding RX J0513. 9- 6 951, which is the obj ect of 
this paper dCowlev et al.lll993l : ISouthwell et al.lll99a here- 
after S96), and for convenience will be hereafter referred to 
as RXJ0513. The SSS exhibit strong emission lines of H, 
Hell and higher ionization species, which probably arise in 
a hot accretion disk or an associated wind. 

Optical spectroscopic and photorn etric studies (e.g. 
iGallanan etlH Il989l: ISmale et alJ 1198811 yielded some bi- 
nary information, mainly orbital periods (in the range ~ 
10 — 24 hrs), yet the nature of the compact object was not 
clear. "Normal" LMXBs containing neutron star or black 
hole accretors radiate in a wider energy range at the typical 
luminosities observed in the SSS. Using the observed SSS 
temperatures and bolometric luminosities, van den Heuvel 
et al. (1992, hereafter vdH92) therefore proposed that these 
sources in fact contain white dwarf compact objects. The 
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Figure 1. Optical light curve of RX J0513. 9-6951. The observations were obtained at SAAO (triangles), CTIO (crosses) and Barfold 
Observatory (diamonds). The times of the start of each X-ray observation are marked. 



near-Eddington luminosities are achieved because the white 
dwarf accretor is able to sustain steady nuclear burning at 
its surface. This requires anomalously high accretion rates 
(when compared, for example, to most cataclysmic variables; 
CVs) of > 10~^Mq yr~^. However, such rates are attainable 
if the donor star is comparable to or more massive than the 
whit e dwarf, so that t hermally unstable mass transfer occurs 
(e.g. iPaczvriskillToTl . 

The vdH92 accreting white dwarf model for the SSS 
still lacks convincing observational confirmation. However, 
circumstantial evidence in its favour was presented by S96 
who showed that the observed bipolar outflow emission lines 
JPakuU et alll993l:ICowlev et alJl99i) in the transient SSS, 
RXJ0513, have velocities comparable to the escape velocity 
of a massive white dwarf. 

RXJ0513 was discov ered in the ROSAT All Sky Sur- 
vey iSchaeidt et aljFlQQSl) . and the extended monitoring of 
the LMC showed its unique characteristic, namely a very 
large variability which made it a type of transient. The ob- 
servations showed that the source brightne ss increased by 
a fact or of ~ 20 during a period of 10 days. ISchaeidt et alJ 
il993l) found that the PSPC spectrum was well-fitted by a 
blackbody of temperature 40 eV, with a column density of 
Nh = 0.94 X lO^'^ cm~^ and bolometric X-ray luminosity of 
~ 2 X 10^* ergs s"\ 

A key breakthrough in our understandin g of the SSS 
came about through the MACHO project iAlcock et alJ 
Il995l) survey of a field in the LMC that contained RXJ0513. 
This led to the production by S96 of an exceptional 3.5 
year lightcurve of the V ^ 17 optical counterpart which re- 
vealed recurrent low states (dropping by ~ 1 mag) at quasi- 
regular intervals (every 100 — 200 days) and remaining low 
for ~ 20-40 days. More importantly, the ROSATX.-ra.y de- 
tections are only reported during such (optical) low states, 
while no outbursts have been observed during the extensive 
optical high states. This showed that the X-ray and optical 
states are anfi-correlated. Furthermore, analysis of the high 



state optical data revealed a low amplitude modulation at a 
period of 18.3 hrs which S96 interpret as the orbital period 
of RXJ0513. 

The remarkable behaviour observed in RXJ0513 is very 
difficult to reconcile with a regular thermonuclear flash 
model which is normally accompanied by radius expan- 
sion and an increased optical luminosity. Instead S96 pro- 
posed a contraction of the white dwarf from an expanded, 
Eddington-limited state to a steady shell-burning phase as 
a cause for the X-ray outburst. This is consistent with all 
the available data, and an X-ray phase due to contraction 
of the envelope during shell-burning has be en es tablished 
observationally for GQ Mus ll6gelma.n et alJll993l) . 

But what is the cause of such contractions? S96 pro- 
posed that this could be due to a decrease in the otherwise 
very high accretion rate, which would also account for the 
simultaneous optical low states. This then draws an analogy 
for RXJ0513 with the VY Scl class of CVs in that they ex- 
perience downward transitions from an otherwise normal on 
state (Honcycutt & Kafka 2004). Of course, the shell lumi- 
nosity will be even higher in the optically bright state, but 
the system is then close to th e top of the s teady burning 
strip in the M — Mwd plane llNomotdll982h which means 
the wh ite dwarf is inflated (perhaps by a factor 3 in radius; 
iKovetz fc Prialnik 199jD. Hence most of the shell luminos- 
ity would then be emitted in the UV or EUV. This leads 
to the prediction that, if this model is correct, as the X- 
ray outburst evolves it will end with its peak emission (or 
temperature) shifting into the EUV and UV. 

We present here the first XMM-Newton observations 
of RXJ0513 through one of its optical low states. XMM- 
Newton is ideal for this study because, via its onboard Opti- 
cal Monitor, we automatically obtain simultaneous UV pho- 
tometry. We use these data to examine the validity of the 
vdH92 paradigm for SSS via the prediction of S96 for the 
X-ray turn-on and turn-off. If we can confirm directly the 
vdH92 paradigm for SSS then this will also establish them 



XMM Spectroscopy of RX J 05 13. 9-6951 3 



as one of the strong c andidates fo r the progenitors of type 
la supernovae (see e.g. lLiviolll996l) . 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Optical Monitoring 

We were granted ToO status for observations of RXJ0513 as 
part of the cycle 2 of XMM-Newton. When the source started 
its ~ 1 mag drop in the optical we would trigger our XMM- 
Newton observations. To enable this we set-up a ground 
based optical monitoring campaign. From 2002 December 5 
(MJD 24520613.5) to 2004 June 11 (MJD 2453167.5) we ob- 
tained V-band and unfiltered observations of RXJ0513 from 
SAAO (Sutherland, South Africa) , CTIO (La Serena, Chile), 
and Barfold Observatory (Victoria, Australia) . Table^gives 
the observing log. Aperture photometry was performed on 
all of the observations using custom packages. The magni- 
tude of the source was determined usi ng the magnitude o f 
a comparison star (V = 16.55; star 2: ICowlev et aljri993fl . 
The optical light curve of RXJ0513 is shown in Figure Q 

On MJD 2452692 RXJ0513 faded by 0.6 - 0.9 mag, 
but by the following day it had re-brightened to its "steady 
level". A few days later it experienced another drop in 
brightness, ~ 0.25 mag, which lasted for ~ 40 days. This 
less pronounced reduction in brightness can be seen in ear- 
lier lightcurves of RXJ0513 (S96), and since the source did 
not dim by > 0.5 mag we did not trigger our XMM-Newton 
over-ride. A characteristic 1 mag drop in brightness did sub- 
sequently occur around MJD 2452980, however we were un- 
able to trigger our observations due to operational restric- 
tions at the XMM-Newton Science Operations Centre pre- 
vailing at that time. Fortunately, RXJ0513 faded again by 
> 0.5 mag on 2004 April 22, and after confirming this drop 
over several subsequent nights, we were able to initiate our 
ToO on April 26. The first XMM-Newton observation oc- 
curred on 2004 April 28, followed by eight more observa- 
tions, with the last being taken on 2004 May 28. The X-ray 
observation log is given in Table |21 The times of the start of 
each X-ray observation are marked in Figure 



2.2 XMM-Newton Observations 



The OM ijMason et alJl2001^ was operated in Rudi-5 imag- 
ing mode during our observations using the UVW2 filter 
(180 — 225 nm). We processed the OM images using the 
omichain routine under the XMM-Newton Science Analysis 
System (SAS) v6.0.0. The photometry was performed with 
omsource using a 5" radius aperture. Background subtrac- 
tion was done using counts extracted from an aperture of 
10" offset from the source, and deadtime and coincidence 
loss corrections were a pplied to the count r ates. 

The EPIC-MO S l|Turner et alJ 1200 1^ and EPIC-PN 
l|Striider et alJl200j ) instruments were operated in timing 
mode, with the thin filter. We reduced the EPIC data with 
SAS v6.1.0. The data were filtered to exclude times of high 
background. The first observation (ID 0151410101) was af- 
fected by high radiation and yielded no useful data. The 
remaining observations were filtered on pattern and energy. 
Since the data were taken in timing mode we filtered the 
data to include only single photon events for the MOS, and 



only single and double photon events for the PN. Only pho- 
tons with energies in the range 0.3 — 10 keV were included. 
The data have been filtered to exclude events that may have 
incorrect energies, for example those next to the edges of the 
CCDs and next to bad pixels. We note that there are no PN 
data for the second observation (ID 0151412101). The ex- 
posure times given in Table |5| reflect the filtering that was 
performed and provide the duration of good data. 

For each EPIC instrument a spectrum for the source 
was extracted over the RAW-X direction, and a region of 
similar size offset from the source position was used to ex- 
tract a corresponding background spectrum. We used an 
extraction region of > 20 pixels which should encompass 
100% of the flux, however, this is uncalibrated. In each case 
we created a photon redistribution matrix (RMF) and an- 
cillary region flle (ARF). The spectra were regrouped by re- 
quiring at least 30 counts per spectral bin. The subsequent 
spectral fitting and analysis was performed using XSPEC, 
version 11.3.1. We extracted source and background light 
curves from the filtered event files, which had been barycen- 
trically corrected, using the same regions. 

The RGS (den Herder et al. 2001) spectra were ex- 
tracted with rgsproc under SAS v6.1.0. Background sub- 
traction was performed with the SAS using regions adjacent 
to that containing the source. Since SAS v6.1.0 automati- 
cally corrects the response matrices for effective area differ- 
ences between RGSl and RGS2, it was possible to combine 
the spectra, including both first and second order data, chan- 
nel by cha nnel for each obs ervation using the method de- 
scribed bv lPaee et alJ l)200,^^ . The resulting combined RGS 
spectrum for each observation was analysed in SPEX 2.00 
tKaastra et al.ii200a') . 



3 UV AND OPTICAL LIGHTCURVES 

We show in FigureElthe OM UVW2 lightcurve of RXJ0513 
with the contemporaneous optical data points obtained by 
Barfold Observatory. The UV intensity is roughly constant 
at the start of the observations (MJD 2453127 - 2453135). 
By MJD 2453137 the UV flux had started to increase, this 
brightening continued during the remaining observations. 
By comparing the start of the rise in the UV intensity to the 
rise in the optical we can see that the UV leads the optical. 
The lightcurves indicate that the UV increase in brightness 
is more gradual than that in the optical. 

We have also plotted the EPIC PN flux in the 0.2-10.0 
keV band (see Section 1^. The overall trend for the X-rays 
during our observing period is a reduction in the X-ray flux. 
There is however a drop in X-ray flux at MJD 2453137 which 
coincides with the start of the rise in the UV flux. The X-ray 
fiux increases after this date, but then continues decreasing. 
The fourth panel in Figure|5|shows the hardness ratio, which 
is seen to decrease as the UV and optical intensity increases, 
see Section |1] for more details. 



4 X-RAY LIGHTCURVES AND HARDNESS 
RATIOS 

Initial inspection of the X-ray spectra indicates that there is 
little emission at energies > 0.8 keV. We therefore extracted 
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Table 1. Observation log of the optical measurements of RX J0513. 9-6951. 



Date MJD Filter Exp Telescope Date MJD Filter Exp Telescope 

+2450000 (s) +2450000 (s) 



2002- 


■12- 


■05 


2613, 


.7858 


V 


30 


SAAO 1.9 m 


2003- 


■03-17 


2715, 


.9919 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■07 


2615, 


,0842 


V 


30 


SAAO 1.9 m 


2003- 


■03-18 


2716, 


.9884 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■08 


2616, 


,0659 


V 


30 


SAAO 1.9 m 


2003- 


■03-19 


2717, 


.9898 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■08 


2616, 


,7776 


V 


60 


SAAO 1.9 m 


2003- 


■03-21 


2719, 


.0053 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■12 


2620, 


,8108 


V 


30 


SAAO 1.9 m 


2003- 


■03-22 


2720, 


.5384 


none 


1800 


Barfold 


14 in 


2002- 


■12- 


■12 


2620, 


,8119 


V 


60 


SAAO 1.9 m 


2003- 


■03-23 


2721, 


.9937 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■13 


2621, 


,9557 


V 


30 


SAAO 1.9 m 


2003- 


■03-24 


2722, 


.9926 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■13 


2621, 


,9563 


V 


60 


SAAO 1.9 m 


2003- 


■03-25 


2723, 


.9901 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■14 


2622, 


,9846 


V 


30 


SAAO 1.9 m 


2003- 


■03-26 


2724, 


.5630 


none 


1260 


Barfold 


14 in 


2002- 


■12- 


■14 


2622, 


.9856 


V 


30 


SAAO 1.9 m 


2003- 


■03-26 


2724 


.9917 


V 


180 


CTIO 1 


.3 m 


2002- 


■12- 


■16 


2624, 


.0331 


V 


30 


SAAO 1.9 m 


2003- 


■03-29 


2727, 


.4295 


none 


1800 


Barfold 


14 in 


2002- 


■12- 


■16 


2624, 


.0343 


V 


30 


SAAO 1.9 m 


2003- 


■03-27 


2725, 


.9853 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■04 


2674, 


.0530 


V 


180 


CTIO 1.3 m 


2003- 


■03-29 


2727, 


.9840 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■05 


2675, 


.0242 


V 


180 


CTIO 1.3 m 


2003- 


■03-30 


2728, 


.9740 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■06 


2676, 


.0441 


V 


180 


CTIO 1.3 m 


2003- 


■03-31 


2729, 


.9751 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■06 


2676, 


.0674 


V 


180 


CTIO 1.3 m 


2003- 


•04-01 


2730, 


.9747 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■07 


2677, 


.0340 


V 


180 


CTIO 1.3 m 


2003- 


■04-02 


2731, 


.9723 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■07 


2677, 


.0832 


V 


180 


CTIO 1.3 m 


2003- 


■04^03 


2732, 


.4571 


none 


1710 


Barfold 


14 in 


2003- 


■02- 


■08 


2678, 


.0512 


V 


180 


CTIO 1.3 m 


2003- 


■04-03 


2732 


.9890 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■09 


2679, 


.0976 


V 


180 


CTIO 1.3 m 


2003- 


■04-04 


2733 


.9926 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■10 


2680, 


.0597 


V 


180 


CTIO 1.3 m 


2003- 


■04-05 


2734 


.9817 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■11 


2681, 


.0393 


V 


180 


CTIO 1.3 m 


2003- 


■04-06 


2735 


.9820 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■11 


2681, 


.5979 


none 


1800 


Barfold 14 in 


2003- 


■04-10 


2739, 


.0329 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■12 


2682, 


.0321 


V 


180 


CTIO 1.3 m 


2003- 


■04-10 


2739, 


.4765 


none 


1440 


Barfold 


14 in 


2003- 


■02- 


■13 


2683, 


.0222 


V 


180 


CTIO 1.3 m 


2003- 


■04-13 


2742, 


.9679 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■14 


2684, 


.0219 


V 


180 


CTIO 1.3 m 


2003- 


■04^15 


2744, 


.4947 


none 


2340 


Barfold 


14 in 


2003- 


■02- 


■14 


2684, 


.0281 


V 


180 


CTIO 1.3 m 


2003- 


■04^19 


2748, 


.9734 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■16 


2686, 


.0188 


V 


180 


CTIO 1.3 m 


2003- 


■04-19 


2748, 


.4462 


none 


1710 


Barfold 


14 in 


2003- 


■02- 


■17 


2687, 


.0335 


V 


180 


CTIO 1.3 m 


2003- 


■04^21 


2750, 


.9637 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■18 


2688, 


.0150 


V 


180 


CTIO 1.3 m 


2003- 


■04-22 


2751, 


.4350 


none 


1800 


Barfold 


14 in 


2003- 


■02- 


■19 


2689, 


.0135 


V 


180 


CTIO 1.3 m 


2003- 


■04-24 


2753, 


.9578 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■20 


2690, 


.0238 


V 


180 


CTIO 1.3 m 


2003- 


■04-26 


2755, 


.4181 


none 


1620 


Barfold 


14 in 


2003- 


■02- 


■21 


2691, 


.0178 


V 


180 


CTIO 1.3 m 


2003- 


■04-27 


2756, 


.9564 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■22 


2692, 


.0747 


V 


180 


CTIO 1.3 m 


2003- 


•04-29 


2758, 


.4204 


none 


1800 


Barfold 


14 in 


2003- 


■02- 


■22 


2692, 


.6165 


none 


1800 


Barfold 14 in 


2003- 


■05-01 


2760, 


.9718 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■23 


2693, 


.0175 


V 


180 


CTIO 1.3 m 


2003- 


•05-02 


2761, 


.9844 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■24 


2694, 


.0247 


V 


180 


CTIO 1.3 m 


2003- 


•05-05 


2764, 


.4242 


none 


1800 


Barfold 


14 in 


2003- 


■02- 


■25 


2695, 


.0212 


V 


180 


CTIO 1.3 m 


2003- 


•05-03 


2762, 


.9541 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■25 


2695, 


.4986 


none 


1800 


Barfold 14 in 


2003- 


•05-07 


2766 


.9501 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■26 


2696, 


.0158 


V 


180 


CTIO 1.3 m 


2003- 


■05-10 


2769, 


.4898 


none 


1440 


Barfold 


14 in 


2003- 


■02- 


■27 


2697, 


.0121 


V 


180 


CTIO 1.3 m 


2003- 


■05-13 


2772, 


.9565 


V 


180 


CTIO 1 


.3 m 


2003- 


■02- 


■28 


2698, 


.0524 


V 


180 


CTIO 1.3 m 


2003- 


■05-14 


2773, 


.9495 


V 


180 


CTIO 1 


.3 m 


2003- 


■03- 


■01 


2699, 


.0354 


V 


180 


CTIO 1.3 m 


2003- 


■05-15 


2774, 


.9516 


V 


180 


CTIO 1 


.3 m 


2003- 


■03- 


■01 


2699, 


.3983 


V 


180 


CTIO 1.3 m 


2003- 


■05-16 


2775, 


.9548 


V 


180 


CTIO 1 


.3 m 


2003- 


■03- 


■03 


2701, 


.0018 


V 


180 


CTIO 1.3 m 


2003- 


■05-17 


2776, 


.9572 


V 


180 


CTIO 1 


.3 m 


2003- 


■03- 


■04 


2702, 


.0000 


V 


180 


CTIO 1.3 m 


2003- 


■05-18 


2777, 


.9529 


V 


180 


CTIO 1 


.3 m 


2003- 


■03- 


■04 


2702, 


.9992 


V 


180 


CTIO 1.3 m 


2003- 


■05-21 


2780, 


.4894 


none 


1980 


Barfold 


14 in 


2003- 


■03- 


■05 


2703, 


.5417 


none 


1890 


Barfold 14 in 


2003- 


■05-23 


2782, 


.4390 


none 


1620 


Barfold 


14 in 


2003- 


■03- 


■06 


2704, 


.0338 


V 


180 


CTIO 1.3 m 


2003- 


■05-23 


2782, 


.9485 


V 


180 


CTIO 1 


.3 m 


2003- 


■03- 


■07 


2705, 


.0100 


V 


180 


CTIO 1.3 m 


2003- 


■05-25 


2784, 


.9439 


V 


180 


CTIO 1 


.3 m 


2003- 


■03- 


■08 


2706, 


.0070 


V 


180 


CTIO 1.3 m 


2003- 


■05-27 


2786, 


.5102 


none 


1800 


Barfold 


14 in 


2003- 


■03- 


■09 


2707, 


.0093 


V 


180 


CTIO 1.3 m 


2003- 


■05-30 


2789, 


.4389 


none 


1800 


Barfold 


14 in 


2003- 


■03- 


■10 


2708, 


.0049 


V 


180 


CTIO 1.3 m 


2003- 


■10-18 


2930, 


.5964 


none 


789 


Barfold 


14 in 


2003- 


■03- 


■11 


2709, 


.0045 


V 


180 


CTIO 1.3 m 


2003- 


■10-26 


2938, 


.6884 


none 


900 


Barfold 


14 in 


2003- 


■03- 


■11 


2709, 


.4714 


none 


1890 


Barfold 14 in 


2003- 


■11-03 


2946, 


.6754 


none 


900 


Barfold 


14 in 


2003- 


■03- 


■12 


2710, 


.0174 


V 


180 


CTIO 1.3 m 


2003- 


■11-08 


2951, 


.6811 


none 


450 


Barfold 


14 in 


2003- 


■03- 


■13 


2711, 


.0027 


V 


180 


CTIO 1.3 m 


2003- 


■12-01 


2974 


.6263 


none 


1800 


Barfold 


14 in 


2003- 


■03- 


■14 


2712, 


.0351 


V 


180 


CTIO 1.3 m 


2003- 


■12-06 


2979, 


.5504 


none 


1710 


Barfold 


14 in 


2003- 


■03- 


■14 


2712, 


.9966 


V 


180 


CTIO 1.3 m 


2003- 


■12-11 


2984, 


.7006 


none 


135 


Barfold 


14 in 


2003- 


■03- 


■16 


2714, 


.9882 


V 


180 


CTIO 1.3 m 


2003- 


■12-13 


2986, 


.6378 


none 


900 


Barfold 


14 in 



The SAAO observations were taken using the 1.9 m telescope with the WRTl CCD. The CTIO observations were 
taJsen using the 1.3 m telescope with the ANDICAM-CCD and Fairchild 447 detector. The Barfold observations 
were made using: the 14 inch @f/6.0 telescoDe with the Audine 401E instrument. 
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Table 1 - continued 



Date MJD Filter Exp Telescope Date MJD Filter Exp Telescope 

+2450000 (s) +2450000 (s) 



ZUUo- 


1 9 


ID 


9QSQ 

zyoy 


Oolo 


none 




JjiiriOlQ 


1 A 


in 


zuu^ 


0*^ 


O'X 
ZO 


'iOS7 
oUo 1 


1 01 


none 


Qfin 
yuu 


Barfold 


1 A 


in 


2003- 


12- 


23 


2996 


4976 


none 


1350 


Barfold 


14 


in 


2004 


04 


16 


3111 


4529 


none 


900 


Barfold 


14 


in 


2003- 


12- 


26 


2999 


3870 


none 


1350 


Barfold 


14 


in 


2004 


04 


21 


3116 


6381 


none 


1350 


Barfold 


14 


in 


2004- 


01- 


08 


3012 


6388 


none 


900 


Barfold 


14 


in 


2004 


04 


22 


3117 


4213 


none 


900 


Barfold 


14 


in 


2004- 


01- 


14 


3018 


6219 


none 


900 


Barfold 


14 


in 


2004 


04 


28 


3123 


4436 


none 


900 


Barfold 


14 


in 


2004- 


01- 


17 


3021 


6854 


none 


810 


Barfold 


14 


in 


2004 


05 


01 


3126 


4854 


none 


900 


Barfold 


14 


in 


2004- 


01- 


28 


3032 


6385 


none 


1350 


Barfold 


14 


in 


2004 


05 


08 


3133 


4803 


none 


900 


Barfold 


14 


in 


2004- 


01- 


31 


3035 


6304 


none 


900 


Barfold 


14 


in 


2004 


05 


12 


3137 


3798 


none 


900 


Barfold 


14 


in 


2004- 


02- 


04 


3039 


6714 


none 


900 


Barfold 


14 


in 


2004 


05 


18 


3143 


5267 


none 


540 


Barfold 


14 


in 


2004- 


02- 


12 


3047 


5078 


none 


720 


Barfold 


14 


in 


2004 


05 


22 


3147 


4069 


none 


900 


Barfold 


14 


in 


2004- 


02- 


19 


3054 


5051 


none 


900 


Barfold 


14 


in 


2004 


05 


29 


3154 


4243 


none 


900 


Barfold 


14 


in 


2004- 


02- 


24 


3059 


5149 


none 


990 


Barfold 


14 


in 


2004 


05 


31 


3156 


4304 


none 


540 


Barfold 


14 


in 


2004- 


03- 


02 


3066 


4980 


none 


900 


Barfold 


14 


in 


2004 


06 


08 


3164 


4391 


none 


900 


Barfold 


14 


in 


2004- 


03- 


14 


3078 


4583 


none 


900 


Barfold 


14 


in 


2004 


06 


11 


3167 


5289 


none 


630 


Barfold 


14 


in 


2004- 


03- 


18 


3082 


5838 


none 


900 


Barfold 


14 


in 
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Figure 2. First panel: OM UVW2 lightcurve of RX J0513.9- 
6951, second panel: the contemporaneous optical data points from 
Barfold Observatory, third panel: EPIC FN lightcurve in the 0.2 — 
10.0 keV band (see Section l5l . fourth panel: hardness ratio (sec 
Section 



source and background light curves from the PN data in the 
energy ranges 0.3 — 0.8, 0.3—0.5 and 0.5 — 0.8 keV. We binned 
the hghtcurves into 300 s bins. We define the hardness ratio 
for our data as (0.5-0.8)/(0. 3-0.5) keV. The light curves in 



the different energy bands and the hardness ratios for each 
observation are shown in Figure 3. 

The X-ray flux is variable over the course of the obser- 
vations. Prominent X-ray dips occur in several observations. 
As the UV and optical brightness start to recover (> MJD 
2453150) the X-ray flux and hardness ratio decrease (see 
Figure |2| fourth panel). The dynamical range of the vari- 
ations in the hardness ratio is smaller than the long-term 
evolution of the hardness ratio e.g. in the second sequence 
of data in Figure 3 (Observation ID 0151412301) two dips 
occur separated by ~ 1.5 hr. 



5 LOW RESOLUTION X-RAY 
SPECTROSCOPY 

For the spectral analysis we have only used the higher signal- 
to-noise PN data. Since no substantial emission is detected 
at energies > 0.8 keV we extracted the PN spectrum over 
0.3 — 0.8 keV for each observation. In each case we fltted 
the PN spectra with a blackbody model modifled by neutral 
photoelectric absorption. 

The initial flts to each spectrum where we allowed the 
equivalent hydrogen column density and the temperature 
to be free resulted in values for Nh and kT which were 
poorly constrained, as expected for such low kT values. 
We modelled all the PN spectra simultaneously and found 
Nh = (0.62+E;;°?) X 10^1 cm'^ and kT = 43.951^^^ eV. 
This value for the equivalent hydrogen column density is 
consistent with the average Galactic absorption in this di- 



rection^ (Blackburn IQq I), Nh = (0.69±0 . 11) x 10^^ cm"", 
and the value inferred bv lGaensicke et al] l(l99^ from HST 
data, Nh = (0.557b0.10) x 10^^ cm"^. We therefore repeated 
the fits to the individual PN spectra, fixing the column den- 
sity at the value found from fitting all of the PN spectra 
simultaneously. 

We find that the blackbody model does not fit the spec- 
tra well (typically xt > 1-9), implying that a more sophisti- 
cated model is required to describe the X-ray emission from 



http: / /heasarc. gsfc.nasa.gov/ftools/ 
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Table 2. Observation log of tiie X-ray measurements of flX J0513. 9-6951. 



Observation 


Date 




MJD 
-1-2450000 


Instrument 


Exp 
(s) 




2004- 


■04- 


■28 


3123, 


.7452 


MOSl 


Q 




2004- 


-04- 


■28 


3123, 


,7452 




Q 




2004- 


-04- 


■28 


3123, 


,7579 


PN 


Q 


01 c;i A^ 91 01 


2004- 


■05- 


■uz 


3127, 


,0789 


MOSl 






2004- 


■05- 


■02 


3127, 


,0789 




10562 




2004- 


■05- 


■uz 


3127, 


,0784 




lUOOO 




2004- 


■05- 


■02 


3127, 


.0784 


RGS2 


16880 




2004- 


■05- 


■uz 


3127, 


,0815 


WiVX C V vv z 


979Q 

Z J Zi^ 




2004- 


■05- 


■02 


3127, 


.1468 


OM TTVW? 


2729 


01 f^l A^ 9901 


2004- 


■05- 


■Uu 


3130, 


.8146 


LVIKJOI- 


1 "^SR^i 

xoooo 




2004- 


■05- 


■05 


3130, 


.8146 


MOS2 


14237 




2004- 


■05- 


■uo 


3130, 


.8274 


PN 


1 Anns 

14:UUO 




2004- 


■05- 


■05 


3130, 


.8141 


RGSl 


17885 




2004- 


■05- 


■uo 


3130, 


.8141 


XV'^OZ 


1 I oou 




2004- 


■05- 


■05 


3130, 


.8173 


OM flVW? 


5000 




2004- 


■05- 


■Uu 


3130, 


.9088 


WiVi C V Vv z 




01 "^141 9^01 


2004- 


■05- 


■09 


3134, 


.9117 


MOSl 


16720 




2004- 


■05- 


HQ 
■uy 


3134, 


.9117 


IVXV_/OZ 






2004- 


■05- 


■09 


3134, 


.9453 


PN 


14049 




2004- 


■05- 


■09 


3134, 


.9111 


RGSl 


17785 




2004- 


■05- 


■09 


3134, 


.9112 


RGS2 


17780 




2004- 


■05- 


■09 


3134, 


.9144 




5000 




2004- 


■05- 


■10 


3135, 


.0267 


OM JIVW? 


3397 


01 "^141 94.01 


2004- 


■05- 


■12 


3137, 


.6860 


MOSl 


6231 




2004- 


■05- 


■12 


3137, 


.6860 


iVJ. v_y O Z 


6443 




2004- 


■05- 


■12 


3137, 


.6987 


PN 


5378 




2004- 


■05- 


■12 


3137, 


.6855 


RGSl 


9^fiS'^ 
zouou 




2004- 


■05- 


■12 


3137, 


.6855 


RGS2 


9^fiS0 
zuuou 




2004- 


■05- 


■12 


3137, 


.6887 


OM JJVW1 


3888 




2004- 


■05- 


■12 


3137, 


.7673 


OM TFVW? 


3397 




2004- 


■05- 


■12 


3137, 


.8403 


OM JJVW1 


3397 




2004- 


■05- 


■12 


3137, 


.9133 


OM JIVW? 


3397 


01 "^1 41 9f^01 


2004- 


■05- 


■J-U 


3141, 


.0642 


MOSl 


lOZOU 




2004- 


■05- 


■16 


3141, 


.0642 


MOS9 


13287 




2004- 


■05- 


1 f\ 
■±o 


3141, 


.0769 


PN 


12192 




2004- 


■05- 


■16 


3141, 


.0636 


RGSl 


13885 




2004- 


■05- 


■16 


3141, 


.0637 


RGS2 


13880 




2004- 


■05- 


■16 


3141, 


.0669 


OM TTVW? 


5000 




2004- 


■05- 


■16 


3141, 


.1584 


OM TFVW? 


3097 


01 "^141 9fi01 


2004- 


■05- 


■18 


3143, 


.3513 


MOSl 


11882 




2004- 


■05- 


■-LO 


3143, 


.3513 


iVl woz 


12043 




2004- 


■05- 


■18 


3143, 


.3641 


PN 


10958 




2004- 


■05- 


1 s 

■-LO 


3143, 


.3508 








2004- 


■05- 


■18 


3143, 


.3508 


RGS2 


13880 




2004- 


■05- 


■lo 


3143, 


.3540 


WiVi LJ V vv z 


ooxy 




2004- 


■05- 


■18 


3143, 


.4319 


OM TFVW? 


4278 


01 "^1 41 9701 


2004- 


■05- 


9fi 
■zo 


3151, 


.2400 


IVXWOX 


1 71 7"^ 
X / X / o 




2004- 


■05- 


■26 


3151, 


.2401 


M099 


17304 




2004- 


■05- 


■zu 


3151, 


.2528 


PN 






2004- 


■05- 


9fi 
■ZD 


3151, 


.2395 




ioUoO 




2004- 


■05- 


■26 


3151, 


.2395 


RGS2 


18080 




2004- 


■05- 


■26 


3151, 


.2427 


OM UVW2 


2819 




2004- 


■05- 


■26 


3151, 


.3090 


OM UVW2 


3278 




2004- 


■05- 


■26 


3151, 


.3806 


OM UVW2 


3278 


0151412801 


2004- 


■05- 


■28 


3153, 


.2732 


MOSl 


15252 




2004- 


■05- 


■28 


3153, 


.2732 


MOS2 


15278 




2004- 


■05- 


■28 


3153, 


.2860 


PN 


14141 




2004- 


■05- 


■28 


3153, 


.2727 


RGSl 


15885 




2004- 


■05- 


■28 


3153, 


.2727 


RGS2 


15880 




2004- 


■0-5- 


■28 


3153, 


.2759 


OM UVW2 


3094 




2004- 


■05- 


■28 


3153, 


.3454 


OM UVW2 


3094 




2004- 


■05- 


■28 


3153, 


.4149 


OM UVW2 


1000 
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Figure 3. The X-ray light curves for the EPIC PN observations of RX J0513.9-6951 in the 0.3 - 0.8 (first panel), 0.3 - 0.5 (second 
panel) and 0.5 — 0.8 kcV bands (third panel). The binning for the light curves is 300 s. The hardness ratio is shown in the fourth panel, 
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"2453130.8274 




0.3 0.5 O.B 

Energy (keV) 

Figure 4. EPIC PN spectra of RX J0513. 9-6951 plotted in 
chronological order of descending time. Each spectrum is labelled 
with its observation time in MJD. At higher energies the spectra 
are dominated by noise and not by real spectral features. 



RXJ0513. At higher energies the spectra are dominated by 
noise and not by real spectral features. We show in Figure^] 
the PN spectra in chronological order. The last two spectra 
show apparent broad absorption features just below 0.5 keV. 
We tried adding one or more absorption edges to the fit (at 
the observed energy they could originate from C V and/or 
C vi); an edge model was unable to reproduce the form of 
the features. The profile of the absorption is too well defined 
to be explained by an edge, at the spectral resolution of the 
PN, and is thus more indicative of discrete line absorption 
(see Section ISJ. 

While the blackbody model is a poor fit to the spectra, 
we can however investigate the evolution of the spectra by 
treating the values obtained as representative of the gross 
spectral properties. In Table^we show the results of the fits 
to the PN spectra. We find that the temperature varies be- 
tween 38.7 — 47.5 eV, and the blackbody luminosity between 
(1.1 — 2.0) X 1 0^^ erg s~^. These va lues are consistent with 
those found bv lSchaeidt et al.l lll993f) from the ROSAT daXa., 
as is the value for Nh ■ We show in Figure |S] the evolution 
of kT, blackbody luminosity, white dwarf radius and X-ray 
luminosity in the 0.2 — 10 keV band, during our observations. 

As the source starts to emerge from the optical low 
state we see that the temperature and X-ray luminosity both 
decrease. This is anti-correlated with the optical and UV 
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Figure 5. The evolution of kT, blackbody luminosity (Lbb), 
white dwarf radius {Rwd) and X-ray luminosity in the 0.2 — 10 
keV band (Lx)- 
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Figure 6. The eight RGS spectra taken during the X-ray high 
state of RXJ0513, in increasing time order from top to bottom. 
The spectra have been offset by 0.25 counts s~^ in the ver- 

tical direction for clarity. 



6 HIGH RESOLUTION X-RAY 
SPECTROSCOPY 

RGS spectra were available for eight observations. All of 
these spectra, which are of very high statistical quality, show 
evidence of deep absorption features and probably some low 
significance narrow emission lines. A plot of all eight spectra 
is given in Figure El 

The absorption lines can be identified with transitions 
from a range of high-ionization ions, including C vi, S xil, 
S XIII, S XIV, Ar XIII, Ar xiv, N vi and N vii, and possibly 
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Table 3. Blackbody fits to the EPIC PN spectra of RX J0513.9-6951. 



MJD 
+2450000 


kT 
eV 


Lbb 

IqSS gj.g g-l 


RwD 
10^ cm 


pa 

lO^^'' erg cm~ s"-"^ 


X 

10~i° erg cm~2 s~l 


10^8 erg s-l 


1038 erg s"i 




3130.8274 


46.281°:?^ 


r, n9+0.08 
^■'-'^-0.02 


1 otr + O.OS 
-1--O3_0.07 


9 qc+O.ll 


C qq + 0.22 

^•'^•^-0.11 


'J- "^-0.03 




6.6 


3134.9453 


47.49t»;g 


^•^^-0.01 






5 02+°-l^ 


69+'' °3 
'^•"^-0.02 


-'-■^'^-0.02 


8.7 


3137.6987 


45 04+''''"' 


1 43+0.08 


1 64+°-"** 


^•'"-0.09 


3.72+0-23 


47+0'04 

-0.03 


1 1^+0.07 
^■^^-0.04 


1.9 


3141.0769 


''S-*' -0.47 


-I yc+0.08 
^•"^-0.02 


1 70+0. 08 
-■-■'^-0.05 




4 57+0.23 




1 q-+0.07 


5.8 


3143.3641 


45 62+"'^* 


1 61+"" 

^•°^-o.a2 


1 69+°-°* 
^■"^-0.06 


-I 0-1+0. 15 


4 20+°-3° 




1 26+°°^ 


4.1 


3151.2528 


39.86l»:«^ 


-I 09+0. 09 
^••^^-0.03 




^•'J^-0.09 


17 + 0.26 

•^■^'-0.13 


09+0.03 

'^•■^^-0.02 


95+° °* 

"■^^"^-0.04 


5.2 


3153.2860 


38.70to;7° 


1 12+"-" 


1 96+°-^2 


n s/i+o-io 

'J-^4_o.07 


9 f,9+0.26 
^■"^-0.14 


'^•^''-0.02 


n 70+0. 08 


3.2 



The value of Nn is fixed at 0.62 X lO^^ cm"^ for all of the fits. Lbb is the blackbody luminosity and RwD is the radius of the white 
dwarf. and are the X-ray fiuxes in the (a) 0.2 — 10.0 keV and (b) 0.1 — 2.4 keV energy ranges, respectively. L'^ and L'^ are the 
X-ray luminosities in the (a) 0.2 — 10.0 keV and (b) 0.1 — 2.4 keV energy ranges, respectively. The errors quoted are at the 90% confidence 
level. 



O VII at 21.6 A. Figure |7| shows three individual features 
plotted in velocity space: C vi Lya at 33.736 A, Ar xill/Ar 
XIV at 27.463 A (blended at the spectral resolution of the 
RGS), and S xiv at 30.441 A. It is clear from this plot that 
the lines are blue-shifted and that there are two or more 
velocity components; the outflow velocity along our line of 
sight is anywhere between zero and ~ 3000 km s~^. 

Figure also indicates that the optical depth of the ab- 
sorbers increases over the course of the eight observations. 
This is seen most clearly in the case of C vi Lya, which is 
barely present in the first spectrum, deepens significantly 
and becomes heavily saturated by the time that the final 
data were taken. The velocity structure of the medium is 
best indicated by the S xiv and Ar xill/Ar xiv features; 
although there are some apparent changes over the course 
of the observations, there appear to be two main velocity 
components at ~ —1000 km s~^ and ~ —3000 km s~^ re- 
spectively. The outflow velocities have not been corrected 
for the systematic vel ocity of the LMC of 280 km s~^ (c.f. 
iHutchines et alj|2002r) . 

Figure |S| shows a comparison of the fifth RGS spectrum 
(MJD 2453141.5636) with a simple photoionized absorption 
model to highlight the positions of the major features. The 
model continuum was chosen to represent the overall shape 
of the first RGS spectrum (MJD 2453127.0784), which has 
the least amount of intrinsic ionized absorption and there- 
fore the most bare continuum visible. The model continuum 
consists of a Compton-scattered blackbody at a temperature 
of 70 eV, with Galactic neutral absorption at the average col- 
umn derived from the fits to the PN data (0.62 x 10^^ cm"^). 
The intrinsic ionized absorption is represented using the xabs 
photoionized absorption model in SPEX 2.00 (Kaastra et 
booa) . which is based on Xstar tKallman fc McCrav.,.,1983) 
output and applies absorption by a column of gas at a given 
ionization parameter ^ and specified (variable) elemental 
abundances. The ionization parameter is defined here as : 



where L is the source luminosity (in erg s~^) in the 
13.6 eV— 13.6 keV range (i.e. the entire photoionizing contin- 
uum), n the gas density (in cm"^) a nd r the source dist ance 
in cm, so ^ has units of erg cm s" jTarter et alJll96gl) . 

The columns of the absorption phases were estimated 
by matching the depth of the modelled C vi Lya feature at 
33.736 A (rest wavelength) with the corresponding feature 
in the data, assuming that carbon is at a Solar abundance 
in the absorbing gas. A gas turbulent velocity of 200 km 
s~^ was chosen (which is below the velocity resolution of 
the RGS; the velocity width of the lines is difficult to es- 
tablish due to the presence of heavily blended features at a 
range of outflow velocities) . The dot-dash line in Figure |S| 
(red in online version) shows the model with all elements at 
Solar abundances. Clearly, the depths of the observed fea- 
tures require far higher relative abundances of at least N, 
S and Ar, and the solid line (blue in online version) shows 
the model with abundances of these elements enhanced so 
that the depths of the relevant features are reproduced. The 
spectrum tails off in the region where oxygen absorption is 
expected, and so it was not possible to gain a useful estimate 
of the relative abundance of oxygen. 

The outflow velocities were set from the observed wave- 
length distributions of the C vi Lya and S xiv (30.441 A) 
features. Ionization parameters of the three phases were then 
estimated by matching the patterns of S absorption in the 
model phases to those in the data (note that the ioniza- 
tion balance of the model was calculated assuming a soft 
excess AGN SED, which is a rather harder spectrum than 
is strictly appropriate for this source). Table 4 gives the 
absorbing columns, ionization parameters and abundances 
relative to Solar of C, N, S and Ar of the final model. We 
emphasize that this model has not been formally fitted to 
the data, and is presented for comparison purposes. 

We note that the columns required to reproduce the 
deep absorption features are one to two orders of magnitude 
higher than the neutral absorption column fitted to the PN 
spectra, which was consistent with the Galactic absorption 
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Figure 7. Absorption features from tlie RGS spectra plotted in velocity space; C VI Ly«, an Ar Xlll/Ar XIV blend and a S XIV line at 
33.736 A, 27.463 A, and 30.441 A, respectively. The dotted lines at km and —3000 km s~^ mark the approximate boundaries of 
the features. 



in the direction of RXJ0513. This is unsurprising since the 
very highly-ionized gas giving rise to the intrinsic absorption 
features requires much higher columns than neutral material 
in order to leave a spectral signature, and, at least in the 
earliest RGS observations, discrete spectral line absorption 
dominates over continuum (edge) absorption. By the time 
the final observations were taken, however, the deep absorp- 
tion features visible even at the PN spectral resolution are 



an indication that the intrinsic ionized absorption column 
has become much deeper. 

It is possible that the apparent highly anomalous abun- 
dances of N, S and Ar are the result of partial covering; if 
the covering factors of the phases are less than unity, but the 
major absorption lines are saturated, increasing the intrin- 
sic column mak es low-abundance ions appear unexpectedly 
prominent (e.g. lArav et alJl2002li . If the relative elemental 
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Table 4. Parameters of the photoionized absorber model compared with RXJ0513 in Figure]^ (sohd line; blue in 
online version). The dot-dash line (online: red) in that Figure has identical parameters except that all elemental 
abundances are Solar. The table lists the following parameters for each absorbing Phase: Nh, the absorbing 
column in units of 10'^^ cm~^; log ^, the log of the ionization parameter where ^ is in erg cm s~^; vtiuc, the 
outflow speed in km s~^; Ac, abundance of carbon as a ratio to its Solar abundance; Ajv, abundance of nitrogen 
as a ratio to its Solar abundance; Ag, abundance of sulphur as a ratio to its Solar abundance; Aat, abundance 
of argon as a ratio to its Solar abundance. 



Phase Nh log 5 Vbiuo Ac An Ns Nat 

1 fO 2.6 800 1 f 2 80 

2 2 2.4 1800 1 20 4 10 

3 4 2.6 2800 1 80 30 90 



abundances in the ionized outflows in RXJ0513 really are 
factors of several above Solar values, this may imply that 
the outflows contain the products of hot CNO burning. As 
stated above, it would be difficult to get an estimate of the 
oxygen abundance from this very soft spectrum due to lack 
of continuum in the relevant spectral region. However, at 
least some of our very rough estimates of the relative abun- 
dances of C and N do seem to be in line with the predic- 
tion (e.g. iLanz et al.ll200^ that carbon should be highly 
depleted with respect to nitrogen in the atmosphere of a 
supersoft source. The prediction that the He/H abundance 
ratio should be enhanced in winds or outflows from white 
dwarfs is not possible to test with this dataset since no spec- 
tral features from these two elements are present in the RGS 
wavelength range, and in any case, they are fully stripped 
at the ionization levels implied by the observed transitions. 

Comparing the models with the observed spectrum in 
Figure |5] it is clear that a simple photoionized absorption 
model is not an adequate representation of the data; this 
is unsurprising since the assumption of absorption from the 
ground state would not be valid in a high density medium 
such as a white dwarf atmosphere. It is also very likely that 
atomic data are missing or perhaps inaccurate in the current 
model for various features in the spectrum; identification of 
ionic species is further complicated by the likelihood that 
the observed features consist of complicated blends of many 
different transitions. Two of the greatest discrepancies are 
in the ranges 25.5 — 27 A and 31.5 — 33 A which correspond 
to the positions of the C vi and C v absorption edges re- 
spectively. 

Approximating the underlying continuum as a black- 
body spectrum is probably still too simplistic; the contin- 
uum and intrinsic absorption will need to be modelled self- 
consiste ntly in order to rea ch a more accurate understanding 
of both. iLanz et al.l ll2005l) have shown that a white dwarf 
atmosphere model is a good representation of the soft X- 
ray spectrum of the supersoft source CAL 83, although they 
found no evidence for a high-speed outflow in that source. A 
white dwarf atmosphere model with the addition of winds 
and outflows may indeed be required to reproduce the RGS 
spectra of RXJ0513. Nevertheless, the present model demon- 
strates clearly that much of the soft X-ray spectral com- 
plexity of RXJ0513 is consistent with the presence of highly 
ionized outflowing gas, and detailed modelling is currently 
underway (Blustin et al., in preparation). 

It is interesting to ask whether the ionized absorber 
model can provide any insights into the nature of the spec- 
tral evolution observed in the soft X-ray band as RXJ0513 



comes out of its optical low state. Further work will be re- 
quired to provide a more certain answer to this, although 
preliminary investigations indicate that much of the varia- 
tion seen in the RGS spectra can be explained by changing 
absorption columns, with the highest columns in the last 
two spectra. The ionization levels, range of outflow veloci- 
ties and elemental abundances appear to be broadly similar 
throughout the series of observations. 



7 DISCUSSION 

The white dwarf photosphere contraction model proposed 
by S96 is based on the prediction that the rise in the X-rays 
occurs after the drop in optical luminosity. From the vdH92 
model, during the optical bright state the accretion rate is 
very high (> 10~'^Mq yr~^). Under these conditions, the 
white dwarf is slightly inflated, and a large majority of the 
shell luminosity is likely emitted in the UV. As the accretion 
rate drops, leading to a reduction in the optical flux, the 
photospher e contracts slightly (e.g. Kovctz & Prialnik 199^; 
|Kat3[l983), this raises the effective temperature and hence 
produces an increase in the X-ray luminosity. 

We do not have X-ray observations before the drop, 
however, previous observations have shown that X-rays are 
not detected when the source is in its optical high state. Our 
observations confirm that while RXJ0513 is in its optical low 
state we see the source in X-rays. 

The behaviour of the UV and optical intensity is 
roughly correlated in the optical low state. However, the UV 
flux starts to rise before the optical luminosity increases. We 
find that as the optical low state progresses the X-ray flux 
decreases, this is anti-correlated with the optical and UV. 
This implies that as the X-ray outburst is evolving the peak 
of the emission has moved into the UV, confirming the pre- 
diction of the contraction model. 

Although our blackbody fits to the PN spectra are poor, 
the values implied by the fits are indicative of the global evo- 
lution of the X-ray emission. We find that the temperature 
and luminosity decrease during the optical low state. The 
radius determined from the fits decreases during the first 
observations. As the optical (and UV) intensity recover the 
radius increases. In the context of the contraction model 
(e.g. Rci nsch et al. ,2000 ; Hachisu fc Katg 2003a.ib|), the en- 
hanced X-ray flux will irradiate the companion star and, 
either by inflating material above the sec ondary's pho to- 
sphere and causing it to be transferred (e.gjKttej|^88|l or 
by heating the magnetic spot region (see lParkeJll979l) . will 
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Figure 8. RGS spectrum number 5 (MJD 2453 141.5636) compared with a simple photoionized absorption model. The dotted line (online 
version: red) shows a version of the model with lAnders & : Grcvcssc' fl98?) Solar abundances, whilst the soUd line (onhne version: blue) 
is a model with enhanced abundances of Ar, N and S with respect to Carbon. The positions of various important features are labelled. 



cause the mass transfer rate to increase again. This coin- 
cides with the white dwarf photosphere re-inflating, hence 
the radius gets larger, and the X-ray flux decreases again. 
Therefore, as RXJ0513 returns to its optical high state the 
source is once again not observable in X-rays. 

Comparing the temperature and luminosity determined 
by the blackbody fits with the white dwarf tracks in the 
Hertzsprung-Russel diagram of Iben (1982; Figure 2) one 
finds agreement with a massive (1.3 — IAMq) white dwarf 
which is somewhat expanded. For a I.SMq white dwarf the 
radius at the stability line is ~ 4.5 x 10* cm. The radius 
inferred from the PN spectral fits i mplies a radius expansion 
by a factor of 3.9. We note that iHachisu fc Katd il2003br) 
derive a larger radius using their model, a factor of ~ 2 
larger than the value we determine from our fits. 

Our spectral fits to the PN data indicate that a more 
complicated model is needed to describe the X-ray emission. 
The higher spectral resolution data from the RGS reveal the 
existence of highly ionized gas outflowing along our line of 
sight at speeds of up to 3000 km s~^, with multiple velocity 
components that evolve over time. The optical depth of the 
ionized absorption appears to increase over the course of 
the eight RGS observations; this is borne out by the lower- 
resolution PN spectra, in which deep absorption features 
become apparent in the final two spectra. The presence of 
this absorption partly explains why the blackbody fits are 
so poor. 

There are a variety of possible causes for the increase in 
optical depth of the intrinsic absorption. We could be view- 
ing the system close to pole-on (i.e. low inclination) through 
the (presumably) bipolar outfiow, and the intrinsic cover- 



ing factor of the outfiow is changing. Possible explanations 
for this include the effects of a magneti c cycle of the opti- 
cal star (c.f. CAL83 and VY Scl stars, lAlcock et al]ll997l: 
iHonevcutt &: K afka"2004. respectively), the expansion of the 
white d warf or the system prec essing with resp ect to our line 
of sigh t JSouthwefl et al.lll997^ . We note that lCowlev et alJ 
i2002l) found a 83.2 d long period in the optical high-state 
MACHO data of RXJ0513 which they attribute to accre- 
tion disk precession. Such a period would be shorter than 
the optical high/low state cycle length. 

Alternatively, it might not be the bipolar outflow we 
are seeing, but actually the process of expansion of the 
white dwarf surface itself, as highly ionized nova by-products 
stream out from the white dwarf surface and gradually build 
into an optically-thick inflated surface, confined by the grav- 
ity of the white dwarf. The increase in optical/UV radiation 
would then occur as the 'intrinsic' X-rays are down-scattered 
in the increasingly optically thick outer infiated layer of the 
white dwarf. The apparently anomalous abundances of Sul- 
phur, Argon and Nitrogen in the ionized outflow may indeed 
imply, if this is not just a covering factor effect due to line 
saturation ( Arav et al. 2002), that the gas originated in the 
nova nucleosynthesis. The abundances may be modified in 
the atmosphere of the steady nuclear burning white dwarf 
via hot CNO burning (c.f. CNO abundances for CAL83 
using NLTE model atmospheres, |Lanz ct al."2005'). From 
the no va evolution model calculations of Priahnk & Kovet3 
il995l) . it follows that a massive white dwarf accreting at a 
high rate of 10~^Mq yr~^ has a recurrence time which can 
be as short as the high/low state cycle length of 100 — 200 
d of RXJ0513. 
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The fact that the observed RGS spectra deviate from 
the blackbody approximation at 24 — 27 A may be due to 
effects of carbon transitions in this range due to a white 
dwarf atmosphere. A white dwarf atmosphere with an effec- 
tive temperature as deduced from the blackbody fits to the 
PN data shows that with increasing gravity (white dwarf 
mass) above the C vi edge (0.49 keV) the atmospheric den- 
sity increases. This shifts the ionization equilibrium toward 
a lower degree of i onization, causing the emission edges to 
turn to absorption iHartmann fc Heiselll996lll997l) . For the 
same gravity the energy of the C absorption (edge) shifts 
to lower energies due to a change of the ionization state of 
C. Thus the C v-vi absorption due to the white dwarf stel- 
lar atmosphere may be at least partly responsible for the 
disagreement found from the blackbod y approximation (see 
iPaerels et al]|200lt iRauch et al1l2005ll . 

The existence of an ionized absorber in the X-ray spec- 
trum of RXJ0513 seems plausible as the absorption lines 
found in the RGS spectrum of the fifth (and other) obser- 
vations of the source require a shift due to a velocity up to 
2800 km s~^ which is consistent with the escape velocity of a 
massive white dwarf. Such absorption lines (but usually with 
lower line-of-sight outflow velocities) are commonly found in 
the soft X-ray spectra of AGN; these were first discovered in 
the Chandra LETGS spec trum of the Seyfert 1 galaxy NGC 
5548 llKaastra et all200(i) and have since been obs erved in a 
range of other sources (see e.g. lBlustin et allbOOSll . The soft 
X-ray absorption lines in the spectra of these AGN originate 
in a low-density photoionized medium. In RXJ0513, how- 
ever, the existence of absorption structures which cannot 
be explained by a photoionized absorber model (especially 
those between 24 — 27 A) indicates that the absorption takes 
place in a very different medium, probably at a much higher 
density. This would be consistent with an origin in a white 
dwarf atmosphere. 

The spectral model we constructed for the fifth RGS 
spectrum in the series does, however, reproduce many spec- 
tral features due to the ionized outflow in RXJ0513; in this 
model, the log ionization parameter ^ of the outflowing gas 
is in the range of 2.4 to 2.6. The ionization structure for 
such an ionization parameter would be consistent with that 
of the outflowing version of model 8 in lKallman fc McCravl 
il982h , which models gas photoionized by a blackbody emit- 
ter. However, it remains difficult to obtain consistency be- 
tween the other parameters inferred by the spectral fits to 
the RGS and EPIC-PN spectra, the hydrogen absorbing col- 
umn density of the ionized gas, the bolometric luminosity of 
the source and the blackbody radius of the source. If we write 
the ionization parameter as ^ = L/NhR where Nh ~ nr^, 
then we obtain with L = 10^^ erg s"\ Nh = 10^^ cm"^ and 
log ^ — 2.6 a value for R = 2.5 x 10^'' cm, which is unrea- 
sonably large for RXJ0513 (the size of the radius inferred 
for the white dwarf in RXJ0513 is 2 x 10® cm). The dis- 
crepancy could be due to the real hydrogen column density 
being much higher than it is in our model. We estimated 
our model absorbing columns assuming that carbon is at its 
Solar abundance in the absorber. If carbon is less abundant 
than this the column density could be underestimated. The 
abundance of carbon in the LMC is only a factor ~ 2 — 3 
less than the Solar abundance which would not have a large 
effect on the column density. The SED used to calculate 
the SPEX xabs absorber model is harder than the apparent 



underlying continuum of this source, and so the ionization 
balance will not be exactl y the same as that of model 8 in 
iKallman fc McCravl (^S^), and the inferred ionization level 
will also be affected if the density is indeed far higher than 
expected in the AGN-type photoionized medium that the 
xabs model assumes (see above). In addition, the outflow- 
ing gas may be ionized by just a fraction of the radiation 
from the nuclear burning white dwarf if it is still close to 
the white dwarf atmosphere. 

We can estimate the mass outflow rate assuming the 
accretion is slightly super-Eddington. The Eddingto n lumi- 
nosity is Lsdd ~ 1.2x10^* M/M q ergs"^ (e. g. Grim m et alJ 
|200i). As Moutf = LEdd/c (see king fc PoundsH2o'o3) . we 
determine with M = IMq, Lsdd = 1-2 x 10^* erg s"^, and 
u = 3000 km s~\ Mont = 2.1 x W'' Mq yr~\ This value for 
the mass outflow rate is consistent with th e wind mass loss 
rate estimated bv lHachisu fc Katd ll2003d) for RXJ0513. 

We can also directly compare the physical parameters 
inferred from the observational data of the low-state egress 
of RXJ0513 with the pa rameters predicted by the model of 
iHachisu fc Katd ll2003bl) . for the high/low state transitions 
of RXJ0513. The evolution of the blackbody temperature in- 
ferred from the EPIC-PN data fol lows roughly the evolutio n 
of the temperature as inferred bv lHachisu fc Katd (l2003lJ) . 
although the temperatures and the radii inferred from the 
observations are about 25% higher and a factor of 2 smaller , 
respectively. One prediction of the lHachisu fc Katd i2003bll 
model are winds from the white dwarf with a wind mass-loss 
rate of ~ 10~^Mq yr~^ during this part of the lightcurve. 
We can estimate the optical depth of such a highly ionized 
wind due to Thomson scattering opacity assuming an out- 
flow velocity of 3000 km s~^ and a radius of 10® cm. For a 
wind mass-loss rate of 10~^A/q yr~^ we calculate an opti- 
cal depth of 1 which is sufficient to obscure the soft X-ray 
source during the low-state egress. 



8 SUMMARY 

We have presented a series of eight XMM-Newton EPIC- 
PN and RGS observations during the late phase (onset of 
egress) of an optical low-state of the supersoft X-ray source 
RXJ0513. Simuhaneous XMM-Newton OM and long-term 
ground based optical monitoring are also reported. We have 
derived the evolution of the fiuxes in the observed bands and 
the evolution of the X-ray spectral parameters from a black- 
body spectral fit to the EPIC-PN data. We find that the 
temperature and luminosity decrease, and an indication of 
an increase in the radius of the blackbody emitter with time. 
During the late phase of the optical low-state we find broad 
spectral dips in the EPIC-PN spectra just below 0.5 keV. 
The RGS spectra show deep absorption features (e.g. C vi, 
Ar xill/xiv, S xiv) which deepen with time and in addition 
some weak emission lines. We model the RGS spectra with a 
Compton-scattered 70 eV blackbody with Galactic neutral 
absorption and additional intrinsic ionized absorption due 
to outflowing gas using a photoionized absorption model. 
Our spectral model for the fifth RGS spectrum, which was 
taken in the deepest part of the optical low state, requires 
velocities of the outfiowing gas of 800, 1800, and 2800 km 
s~^, ionization parameters log ^ of 2.4 to 2.6, and hydro- 
gen columns of at least 2 x 10^^ to lO^'^ cm~^. We find that 
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the spectral model lies significantly above the observed spec- 
trum in the 24—24.5 A and 25.5 — 27 A regime which is likely 
to be due to effects related to the white dwarf atmosphere. 
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